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Abstract 

We present a theoretical study on the spin-dependent transport through a spin valve consisting 
of graphene sandwiched between two magnetic leads with an arbitrary orientation of the lead mag- 
netization. No gate voltage is applied. Using Keldysh's nonequilibrium Green's function method 
we show that, in absence of external magnetic fields, the current- voltage curves are nonlinear. 
Around zero bias the differential conductance versus bias voltage possesses a strong dip. The 
zero-bias anomaly in the tunnel magnetoresistance (TMR) is affected strongly by the leads spin 
polarization. Depending on the value of the bias voltage TMR exhibits a behavior ranging from 
an insulating to a metallic-type. In presence of a static external magnetic field the differential con- 
ductance and TMR as a function of the bias voltage and the strength of the magnetic field show 
periodic oscillations due to Landau-level crossings. We also inspect the effects of the temperature 
and the polarization degrees on the differential conductance and TMR. 

PACS numbers: 75.47.-m,85.75.-d ,81.05.Uw 
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I. INTRODUCTION 



Recent advances in nanoscience techniques opened the way for the creation and the in- 
vestigation of the two-dimensional carbon, also called graphene [H[2|0]0l- This system is 
a monolayer of carbon atoms packed densely into a honeycomb lattice, and can be viewed 
as the basic building block for many carbon-based materials with other dimensionalities, in- 
cluding fullerene, nanotube and graphite, etc. Its low energy band structure consists of two 
inequivalent pairs of cones with apices located at the Brillouin zone corners[5i]. In these cones, 
the energy-dispersion relation is linear, and the dynamics of the charge carriers is governed 
by a massless Dirac-type equation. The form of the electronic band structure is expected 
to lead to a number of unusual electronic properties in graphene such as the anomalous 
quantized Hall effect, the absence of the weak localization and the existence of the mini- 
mal conductivity [6]. Graphene is also an interesting candidate for transport-applications, 
in particular for spintronics:The mobility is remarkably high and the carrier density is con- 
trollable by agate voltage. In addition, spin-dependent interactions can be exploited for 



ll|[l2|. Motivated by these facts, the 



the control of the magnetoconductance [7[[8|j^|l0j 
spin-dependent properties of graphene are in the focus of current research; e.g., E. W. Hill 
et al. fabricated graphene spin-valve device and observed a 10% change in the resistance 



as the electrodes switch from a parallel to an antiparallel state 



13|. Recent experiments on 



spin injection in a single layer graphene show a rather long spin-flip length(^ 1/im) at room 
temperature |U]. Spin injection into a graphene thin film has been successfully demonstrated 
by using non-local magnetoresistance measurements W.H. Wang et al. measured 

the magnetoresistance of mesoscopic graphite spin valve devices and observed a cusp-like 
feature of the magnetoresistance versus the applied bias and pointed out the importance of 
spin-dependent interfacial resistance for spin injectionjl?!. 

In this work we investigate theoretically the spin-dependent transport through a graphene 
spin valve device with ferromagnetic leads having arbitrary spin-polarization directions. No 
gate voltage is applied. Utilizing Keldysh's nonequilibrium Green's function methodjisl we 
calculate the density of state (DOS) and the electrical current in the ferromagnet-graphene- 
feromagnet(FM-G-FM) coupled system. The differential conductance and the tunnel mag- 
netoresistance are also calculated without and with a static external magnetic field at finite 
temperatures. We found that at zero magnetic fields, the current- voltage curves in this spin- 
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tronic structure show a nonlinear characteristic, the differential conductance as a function 
of the applied voltage exhibits a strong dip near zero bias. The behaviour of the zero-bias 
anomaly in TMR is governed by the leads spin polarization. With increasing the tempera- 
ture, the dip in the differential conductance and the cusp in TMR near zero bias diminish. 
When both the spin-polarization degrees and the relative angles of the two ferromagnetic 
moments are large, the differential conductance is small due to the influence of both the DOS 
in graphene and the conventional spin- valve effect. In the presence of a static magnetic field, 
the differential conductance and TMR show periodic oscillations due to a resonant trans- 
port though the Landau levels when the bias-voltage values are appropriate. At zero bias 
voltage, the differential conductance versus the temperature show a behavior different from 
the field-free case. We attribute this fact to the breaking of the insulator-type properties of 
graphene at finite magnetic fields. 

The rest of this paper is organized as follows: In section 2, we introduce the model 
and derive the current formula in the absence of the magnetic field. In section 3, the 
magnetotransport properties of this system are computed at finite external magnetic fields. 
In section 4, the corresponding numerical results are given. Finally, a summary is presented. 

II. THEORETICAL MODEL 

We consider a spin valve device consisting of a graphene layer contacted to ferromagnetic 
electrodes, as shown in FiglH The moment M^, of the left electrode is assumed to define the 
z direction, while the moment of the right electrode deviates from the z direction by a 
relative angle Q. A bias voltage V is applied between the left and the right electrode. The 
electrical current flows in the x direction. The left and the right electrodes can be described 
by Hamiltonians 

i^L = Xl^kLa^i^CkL^, (1) 
k,cr 

Hr = '^[£R(k) - aMn cos 6']Ck^^CkR<x - Mr sin 6*4^^ Cki?^, (2) 

k,cr 

where ekoo- is the single electron energy associated with the momentum k, the spin a and 
the a electrode. c]^a^{ckaa) creates (annihilates) an electron with the energy ekao-- 
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The tight-binding Hamiltonian of the electrons in graphene is given by 

i7G = -t $^(4A- + H.c.), (3) 

where al^{ai^f^) creates (annihilates) an electron with the spin a on the position Rj of the 
sublattice A, h\^{hiu) creates (annihilates) an electron with the spin a on the position 
on the sublattice B, and t is the nearest neighbor hopping energy. In the momentum 

space the Hamiltonian Hq is rewritten as 

Hg = 5^[0(q)a;,6q<, + 0(q)*6i„aq,], (4) 

q,(7 

where 

3 

(j){ci) = -tJ2e'''-'-' with 5i = |(l,v^,0), 52 = |(l,-v^,0), 53 = a(l,0,0). 

i=l 

Here a is the lattice spacing. Diagonalizing the Hamiltonian (jlj) one finds 

E±(k) = ±t|0(k)|, 

which can be linearized around the K points of the Brillouin zone leading to the dispersion 

E±{k) = ±VF\k\, (5) 



where vp = 3ta/2 is the Fermi velocity of electron [t ~ 2.3eV 19|) . The coupling between 
the electrodes and graphene is modeled by 

Ht = ^Y1 [^kaqcLa«q- + H.C.], a = L,R. (6) 



kqcKj 

Tkaq is the coupling matrix between the a electrode and the graphene; is the number of 
sites on the sublattice A. 

The electrical current from the left electrode to the graphene sheet is obtained from the 
time evolution of the occupation number operator of the left electrode, i.e. 

I = e{A) = j{[H,AfL]), Ar, = J]4,^Ck,.. (7) 

kiT 

Using the nonequilibrium Green's function method, Eq.([7]) can be further expressed as 

qq' 



where Tr means the trace in the spin space and is the Fermi distribution function at 

the energy e. 

Gt:Sait-t')=^{a\,At'M)) 
is the matrix expression for the lesser Green's function, ^^^^/^(e) and G'^^^^^{e) are 2x2 
matrices in the spin space in the sublattice A describing respectively the retarded and the 
advanced Green's function. The line width matrix F^qq/ is given by 

/ , \ 

where TkLq stand for the coupling of graphene to the electrodes. To evaluate I from Eq.® 
the retarded Green's function G'[^^,^^{e) needs to be calculated. Here we consider electrons 
near the Fermi level which contribute predominantly to tunneling. In this case one may 
assume the coupling matrix TkLq to be independent of q and set F^qq, = F^. Standard 
Green's function technique [ISi] delivers 

In the above expressions we introduced 



gqLai.!^) = / , ■ NO 1 ,/ M2 > ^ud T{e 



where 



with 



e±^r?)2- 10(g) |2' l-raam^ie) 



&) = ^ E^^«,c,a(^), and J:'-'%e) = T-[rL{e) + RTRm^] 
q 



cos 2 — sm 2 



sm 2 cos 2 



R = 

In Eq.®, G^^^,^{e) can be derived by applying Keldysh equation 

^qa,q'a(^) = ^qa,q'a(^)^^ (^)^qa,q'a(^) 

with 

E<(£) = t[fLie)TLie) + fRie)RTRi6)R^]. 

Substituting the expressions of the graphene Green's functions in Eq.([8]), and after a straight- 
forward calculation we obtain the tunneling current as 

I = t/tTr{[^:(e)(i?F^(e)i?t)]^a(,)r^(,)}[/^(,)_/^(,)], (11) 



where 



GT{e) = J:Y: GTa.'a{e) = , JZ^L... (12) 



Introducing a cutoff kc leads to 

ra:ie) = -Foie)Ttnpois), (13) 

where 

Foie) = ^\n Po(^) = ;^^^p-kl). (14) 

D = Vpkc stands for a high-energy cutoff of the graphene bandwidth. Invoking the Debye's 
prescription, we choose kc such that the total number of states in the Brillouin zone is 
conserved after the linearization of the spectrum around the K point. Hence, Eq. (fT3|) is 
accurate for e <^ D, i.e. e has to be in the region where the linearization of the spectrum is 



justifiable which is roughly estimated |2(| to be [-1.6 eV, 1.6 eV]. In Eq. (fT2l) we assumed a 
symmetrical voltage drop as hl^r = EF±^eV, and put Ep = in the numerical calculations. 
The TMR at the angle 6 is conventionally defined as 

TMRiS) ^ m_m. (15) 

III. MAGNETOTRANSPORT AT FINITE EXTERNAL MAGNETIC FIELDS 

In the presence of a static external magnetic field, the description of the transport prop- 
erties of electrons in a honeycomb lattice becomes much more involved due to the coupling 
between graphene and the electrodes associated with the Hofstadter problem. To circum- 
vent this situation we describe the electrons in the honeycomb lattice as Dirac fermions in 



the continuum. At first we introduce the field operators 



2l||22| 




where Ib = 1/VeB is the cyclotron length, n = 0, 1,2, a = ±1, and (f)n{x) is the nth 
eigenfunction of the usual one-dimensional harmonic oscillator. The Hamiltonian describing 
the electrons in graphene acquires the second-quantized form 

Hg=Y^ E{n, a)dl^c.^dgnaa (17) 

qnaa 



where 

E(n, a) = aUcVn + 1 



is the Landau level with Uc = vpV^eB. The sum over integer n's is cut off at A/". The 
coupling between graphene and ferromagnetic electrodes is 

kqcr 

kqnaa 
kqa 

+ E [(i>n{L - qlV) + a(pn+l{L - qlV)] [TkWjcl^dqnaa + T*^d\^^„CkRa\. 

kqnaa 



(18) 



where L is length of graphene. Similar to the calculation of Eq. ljTTj) we obtain the electrical 
current in the form 

Jl = -n! tTr{[e - (1 + £X)Sr'(l + eXfRTRR^[e - (1 + eX)i:'^Y'T L}[f R{e) - 

(19) 

where 

^ ^ = — ^? — V = — ^H"^ - ^ ^ -2^ 

^-^ e — E(n, a) uJt uji ojf. 

ly^ rM \ ' / L. (_. 



with '^{z) denoting the digamma function 23|. The electric conductance can be obtained 



by dJi/dV . For a small bias voltage we obtain the Landauer-Biittiker-type expression 

G = jT^ff, (20) 

where 

T,ff = Tr{[EF-{l + EFX{EF))T.X\l + EFX{EF)fRTRR^[EF-{l + EFX{EF))T.X^TL} 

plays the role of an effective energy-dependent transmission coefficient. In the limit 5^0, 
we further obtain for 6* = 0, 

^^ff = llEl + \X{Tl + Tir ^^^^ 
where A = {2J\f + 'if + 2{J\f + 1){U + 2){2J\f + ?>)^. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



Adopting the wide bandwidth approximation for the graphene spin valve system we ne- 
glect the energy dependence of the line width functions T'^{e). Denoting the spin polarization 
of the left and the right electrodes by respectively pl and we write T]^^ = T]^^ = ro(l±p), 
here Fq describes the coupling between the graphene and the electrode without the internal 
magnetization. Here we assumed the left and the right electrodes to be of the same material. 
In the following numerical calculation, we take Fq as the energy scale. We calculate the DOS 
in graphene via relation p{e) = — Ga'^'^i^)- Fig l2] shows the DOS as a function of en- 

ergy for different polarization p and magnetization angle 6. It is clearly observed that the 
DOS displays a dip structure with the energy. For nonzero energy, the DOS increases with 
increasing p, however decreases with increasing 6. This is caused by the different tunneling 
rates for up and down spins owing to the s plit ting of DOS of the ferromagnetic leads. This 
splitting acts as an effective magnetic field 2^ reaching values much larger than externally 
applied magnetic field 251]. This results in a spin dependence of the DOS in the central 
region. While at the zero energy point, the DOS is independent of p and 6. This stems from 
the nature of Dirac point in the graphene. 

The bias dependence of the electrical current and the differential conductance G = dl/dV 
are shown in Fig l3] for parallel electrodes magnetizations and for the different polarization p, 
The nonlinear behaviour of G and the strong dip at zero bias are in line with the experimental 

n 

observations p/7| , and are at variance with the typical behaviour when a Fermi liquid were in 
central region instead of graphene. In that case the electrical current is proportional to the 
applied voltage at small bias and obey the Ohmic law 261]. The results obtained here resemble 
rather the ferromagnet-insulator-ferromagnet(FM-I-FM) junctions 27| 28|]. The reason is 
that: the DOS in graphene diminishes at the Fermi level. Hence, graphene sheet can be 
viewed as a tunneling barrier at the zero energy point, similar to FM-I-FM system. With 
increasing the polarization p, the portion of spin-up electron states increases while for the 
spin-down electrons decreases. However, the scattering of the former is larger than that of 
the latter, thus we conclude that the conductance G decreases with increasing p at nonzero 
bias (cf. Fig|3] (b)). While G is almost independent of p at zero bias, which stems from 
the fact that the spin transport through the Dirac point of graphene is ballistic due to its 
insulator-like properties. 
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The bias dependence of the conductance at different temperatures T and angles 9 are 
shown in FigJH The conductance is roughly independent of T at large bias. Near 1^ = 
the conductance increases with increasing T . This behavior is also in contrast to usual 
Fermi liquids where the conductance decreases with increasing temperature because ther- 
mal fluctuations enhance the scattering of conduction electrons and thereby contributes to 
the resistance of system 26|]. In our case, the graphene is equivalent to a barrier at the zero 
energy point. Near zero bias voltage, the thermally excited electrons are dominant in tunnel- 
ing process. Therefore, with increasing temperatures, the increase of the thermally excited 
electrons enhances the conductance. On the other hand for large bias, the contributions 
to the conductance stem mainly from electrons with excess energies well above the Fermi 
level (as dictated by the applied voltage) which leads to a very weak dependence of the con- 
ductance on the temperature at large bias. The conductance as a function of 6 (FigH] (b)) 
follows the conventional behaviour of magnetic junctions such as the ferromagnet-quantum 
dot-ferromagnet system 29|] 301] 311] • When 9 changes from to tt the number of spin-up 
and spin-down electrons is rearranged. Figl5] shows TMR ratio as a function of the applied 
voltage for different polarizations p and temperatures T. A pronounced cusp-like feature 



appears at zero bias in line with experimental observations [17|. We assign this behaviour to 
the result of a non-trivial combined effect of graphene and conventional spin-valve proper- 
ties, evidenced by the dependence of TMR on the polarization p at a fixed bias voltage (cf. 
FiglSl^a)). When increasing the polarization p the contribution of spin- up states relative to 
the spin-down is increased resulting in an increase of TMR for the entire bias range. How- 
ever, the TMR changes in a non-linear manner: the TMR value at zero bias becomes larger 
than that at nonzero bias. This is because for a small bias graphene behaves as an insula- 
tor. In contrast, for a large bias graphene behaves more like a metal, which is an essential 
difference between graphene and other materials in the central region. The spin tunneling 
is ballistic through an insulator in contrast to a metal, whence the TMR is more enhanced 
around zero bias. There is almost no changes with the temperature T for large bias (see 
Figl5](b)). With increasing temperatures more thermally excited electrons contribute to the 
electrical currents for the parallel and the antiparallel configurations, in fact, the increase of 
the electrical current is faster for the antiparallel configuration. Therefore, we can conclude 
that the TMR decreases at higher temperatures and the zero-bias anomaly diminishes in 
this situation. 
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The bias and the magnetic field dependencies of the differential conductance G for the 
different temperature in the parallel configuration of magnetization is shown in Fig|6l G 
versus the applied voltage exhibits an oscillating behaviour. Each conductance peak corre- 
sponds to resonant transport through a Landau level. As bias voltage increases, the distance 
between two neighboring peaks decreases due to the decrease of the distance between neigh- 
boring Landau levels. The differential conductance oscillates as a function of the magnetic 
field, as shown in Figl6] (b). This is due to resonant transport through the Landau levels at 
a particular magnetic field strength and the applied bias values. With increasing tempera- 
tures (Figj6](a)) the differential conductance deceases for all bias voltages. In particular, at 
zero bias the temperature dependence of the differential conductance is different from the 
magnetic field-free case. The explanation for this phenomena is as follows: the magnetic 
field lifts the insulator behavior at Dirac point in graphene, and thus thermal fiuctuations 
suppress the conductance. FigJT] shows the TMR as a function of the bias voltage and the 
magnetic field strength for different polarizations p. TMR reaches a minimal value at bias 
voltages corresponding to the conductance peaks and increases with increasing p which is 
nothing but a spin- valve effect. 



V. SUMMARY 



In conclusion, we have studied the spin-dependent transport through a graphene spin 
valve device for a non-collinear configuration by means of Keldysh's nonequilibrium Green's 
function method. It is found that at a zero magnetic field, the current-voltage curves show a 
nonlinear behaviour. The corresponding differential conductance exhibits a strong dip near 
zero bias. The TMR shows a zero-bias anomaly that depends on the leads-spin polarization. 
Increasing the bias TMR follows a behaviour akin to a metallic or an insulating system 
depending on the value of the bias. In the presence of an external magnetic field, the 
differential conductance and TMR oscillate periodically due to a resonant transport through 
Landau levels. At zero bias the differential conductance versus the temperature reveals a 
behavior different from the magnetic field- free case. 
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FIG. 1: Schematic illustration of the system considered in this work, the graphene is connected to 
two magnetic leads by the tunneling barriers. The moments of the leads are aligned by a relative 
angle 6, and the coupling matrix between a{a = L,R) electrode and graphene is Tka- 
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FIG. 2: DOS as a function of energy e for different polarization p sX = Q {&) and for different 
angle 9 at p = 0.6 (b). The other parameters are taken as-D = 2ro and B = 0. 
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FIG. 3: The bias dependence of the electrical current / (a) and differential conductance G (b) for 
different polarization p aX 9 = Q and ksT = COOSFo. The other parameters are taken the same as 
those of Figlll 
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FIG. 4: The bias dependence of the differential conductance G for different temperature T at = 
(a) and for different angle 9 at ksT = O-OOSFo (b). The parameters are taken the same as those of 
FigEl 
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FIG. 6: The differential conductance as a function of the bias voltage at /xsB = IFq (a) and of the 
magnetic field at eV = SOOFq (b) for different temperature T in the parallel configuration. The 
parameters are taken as p = 0.6 and M = 100. 
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FIG. 7: The TMR as a function of the bias voltage at fj-sB = ITq (a) and of the magnetic field at 
eV = SOOFo (b) for different polarization p at ksT = ClFo. The other parameters are taken the 
same as those of Figl6] 
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